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COMP.VRISON  OF  Wf  SEISMIC  MEASUREMENTS  TAKEN  IN  THE  DIAMOND  DUST, 

DIAMOND  MINE,  AND  MINE  DUST  HE  SERIES  OF  EVENTS 

1 

INTRODUCTION 

This  FOR  complements  the  previously  Issued  technical  report  (ref.  5)  for 
the  MINE  DUST  HE  event  and  contains  the  final  results  of  the  analyses  which  have 
been  performed  on  all  of  the  MIGHTY  MITE  events  but  which  have  not  been  re¬ 
ported  previously  (refs.  3,  4,  5). 

The  Geophysics  Section  of  the  Infrared  and  Optics  Division,  Willow  Run 
Laboratories,  University  of  Michigan  established  six  portable  seismograph 
stations  to  record  this  event.  The  same  close-in  locations  used  for  the  pre¬ 
vious  events  of  this  series  were  reoccupled  (fig.  1).  In  addition,  a  station 
was  established  near  Meadsvllle,  Arizona  about  225  km  southeast  of  ground  zero. 

Stations  of  the  Nevada  seismograph  net  were  also  used  through  the  co¬ 
operation  of  the  National  Oceanic  and  Atmospheric  Administration  (NOAA) ,  and 
the  United  States  Geological  Survey  (USGS).  Figure  2  is  a  map  showing  the 
locations  of  stations  in  the  various  nets. 

2 

FIELD  MEASUREMENT  DATA  AND  RESULTS 

2.1.  UNIVERSITY  OF  MICHIGAN  PARTICIPATION 

The  Geophysics  Section  of  the  Infrared  and  Optics  Division,  Willow  Run 
Laboratories,  University  of  Michigan  established  six  portable  seismograph 
stations  to  obtain  data  for  the  MINE  DUST  HE  event.  The  five  close-in  sta¬ 
tions  (fig.  1)  were  located  in  the  same  positions  as  used  for  DIAMOND  DUST, 
DIAMOND  MINE,  and  the  DIAMOND  MINE  HE  shot  of  February  4,  1971.  All  stations 
were  equipped  with  matched  three- component  short-period  seismometers  and 
magnetic  tape  recorders.  An  additional  high  frequency  vertical  seismometer, 
a  low  frequency  seismometer,  and  a  microphone  were  used  at  the  Pahute  Mesa 
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Figure  1.  Locations  of  seismograph  stations  established 
by  the  University  of  Michigan  for  the  MIGHTY 
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Figure  2.  NOAA/EBL  CGS  Net  i 
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site  located  about  6  km  north-northeast  of  ground  zero.  A  sixth  station 
equipped  with  a  three-component  short— period  seismometer  and  a  tape  recorder 
was  located  near  Meadsville,  Arizona  at  a  distance  of  about  225  km  southeast 
of  ground  zero. 

A  summary  of  liformatlon  about  the  sites  appears  in  Table  1.  The  hori¬ 
zontal  distances  and  bearings  are  with  respect  to  the  MINE  DUST  HE  \7orking 
point  and  magnetic  North.  The  values  were  obtained  from  Holmes  and  Narver 
at  NTS  who  did  the  s’.rveying.  The  values  for  the  North  site  are  only  an  esti¬ 
mate.  The  surveyors  located  a  point  which  was  northwest  of  the  North  site. 

The  elevations  are  estimated  from  Department  of  Interior  Geological  Survey 
maps  (Tipppah  Spring  and  Mine  Mountain).  The  geology  of  the  stations  are  a 
result  of  on-site  inspections,  the  above  mentioned  geological  maps  and  seis¬ 
mic  refraction  profiles  conducted  as  described  below. 


TABLE 

1.  RECORDING 

STATIONS 

Station 

Bearing 

Distance 

(km) 

Elevation 

(ft) 

Station  Geology 

South 

183 ‘’57' 

2.6A 

5190 

Upper  Miocene  or 
Pliocene  Tlva 

Canyon  Member  of 
Paintbrush  Tuff 

Southeast 

119 “59' 

2.59 

5005 

Paintbrush  Tuff 

East 

68“10' 

2.A8 

4890 

Alluvium  over 

Miocene  Tuff,  Indian 
Trail  formation  (older 
than  Paintbrush) 

Northeast 

12“0' 

2.82 

5250 

Pa.  or  Permian 

Tippipah  Limestone 

Pahute  Mesa 

21°30' 

6.05 

4794 

Alluvium 
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Reversed  refraction  profiles  were  made  for  each  site.  Three-component 
seismometers  were  located  at  each  end  of  a  76  m  line  and  the  seismic  source 
was  moved  along  the  line  at  eleven  equally  spaced  locations.  The  line  ’.an 
from  the  site  back  toward  the  working  points.  Thus,  for  example,  the  North 
site  used  for  the  events  was  at  the  north  end  of  the  refraction  profile  for 
that  site,  A  summary  of  the  results  are  given  in  figure  3  and  the  accompany¬ 
ing  table.  They  were  derived  using  standard  refraction  techniques  (Dobrin, 
pp.  81-84),  The  two  northerly  sites  and  the  East  site  have  unconsolidated 
dry  sand  at  the  surface  that  accounts  for  the  low  velocities.  The  other  two 
sites  have  bedrock  at  or  very  near  the  surface  and  at  the  South  site  the  sur¬ 
face  rock  was  badly  fractured.  This,  along  with  the  cracked  regions  from 
previous  shots  (fig.  4),  caused  lower  transmission  velocities.  In  addition 
to  the  layers  revealed  in  these  prof  ies ,  there  is  another  below  them  which 
has  a  velocity  in  the  neighborhood  of  4  km  per  second.  This  was  inferred  from 
the  first  compressional  arrival  time  of  the  various  events.  The  greater  part 
of  the  propagation  path  is  probably  through  the  Indian  Trail  Formation  tuff 
beds. 

Good  data  were  obtained  at  nearly  every  site  for  all  events.  However, 
for  DIAMOND  DUST  the  transverse  component  at  the  East  site  was  s'’tghtly  over¬ 
recorded  and  the  tape  recorder  speed  at  the  North  site  was  not  constant.  The 
East  site  overrecorded  data  are  considered  unreliable  above  10  Hz  and  the 
North  site  data  considered  unreliable  below  10  Hz,  WWVB  was  not  recorded  at 
the  East  site  for  DIAMOND  MINE  HE  and  hance,  could  not  be  used  in  analysis 
where  timing  was  critical.  The  output  from  the  vertical  seismometer  at  the 
Southeast  site  for  DIAMOND  MINE  gave  evidence  of  malfunction  and  was,  therefore, 
not  included  in  the  data  analyzed.  All  other  recordings  were  good. 

The  azimuthal  asymmetry  noted  in  Willis  and  Hand  (1971)  for  the  DIAMOND 
DUST  event  was  observed  for  the  other  three  events  as  well;  the  wave  train 
duration  increases  and  the  frequency  content  decreases  from  north  clockwise 
to  south.  The  East  site  consistently  had  about  double  the  signal  level  which 
would  be  expected  based  on  the  data  gathered  at  the  other  sites.  At  the 
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P-wave  velocity 

km/sec  i 

I  ct2  degrees 


Pahute  Mesa 

7.6 

.363 

1.05 

3.0 

North 

0.9 

.290 

1.22 

-0.5 

East 

15.5 

.472 

1.62 

1 

2.5 

Southeast 

tau 

.671 

1.10 

7.3 

South 

18.0 

.655 

1.15 

5.2 

Figure  3.  Summary  of  refraction  survey  results  at 
the  five  closest  sites  occupied  for'  the 
MIGHTY  MITE  series. 
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Pahute  Mesa  site,  except  for  high  frequencies  for  DIAMOND  MINE,  the  signal 
levels  were  larger  than  expected  on  the  same  basis.  Geological  explanations 
have  been  given  for  these  observations  in  the  above  mentioned  reports. 

Detailed  spectral  analyses  of  the  MINE  DUST  HE  data  were  made  for  the 
five  close-in  stations.  Plots  of  these  spectra  appear  in  Appendix  A.  For 
comparison  purposes  with  the  earlier  three  events  of  this  series  composite 
plots  of  ground  particle  velocity  vs.  frequency  were  made  for  all  four  events, 
Figures  B-1  through  B-10  (Appendix  B)  contain  the  spectral  analyses  of  the 
first  compressional  wave  arrival  recorded  at  each  site  on  the  vertical  and 
longitudinal  component  seismometers.  Particle  velocity  spectral  ratios  of 
MINE  DUST  HE/DIAMOND  MINE  HE  as  a  function  of  site  and  wave  type  are  given 
in  Appendix  C.  Similar  curves  for  DIAMOND  DUST/DIAMOND  MINE  are  shoimin 
Appendix  D.  I  :  can  be  seen  (Appendix  B)  that  the  MINE  DUST  HE  shot  has  sig¬ 
nal  levels  slightly  larger  than  the  DIAMOND  MINE  HE  shot  at  most  frequencies 
but  that  the  relationship  is  frequency  dependent.  In  order  to  quantify  this 
relationship,  particle  velocity  ratios  were  computed  for  the  first  compres¬ 
sional  (P)  wave  arrival  and  the  maximum  shear/surface  wave  (S  )  recorded 
at  each  site.  Geometric  means  for  these  ratios  were  computed  using  that  por¬ 
tion  of  the  spectrum  which  contained  the  predominant  signal  levels  and  which 
were  well  above  background  noise.  These  mean  ratios  are  summarized  in  Table 
2.  The  overall  mean  MINE  DUST  HE  vs.  DIAIWND  MINE  HE  particle  velocity  ratio 
was  1.287.  Also  shown  in  this  table  are  the  means  of  the  vertical  and  longi¬ 
tudinal  components  for  first  P  from  all  stations  and  all  three  components  for 

maximum  S.  The  signal  level  associated  with  the  latter  is  usually  the  largest 
signal  on  the  seismogram. 


TABLE  2. 

AVERAGE  SPECTRAL  RA'i  lOS  MINE  DUST  HE  vs.  DIAMOND  MINE  HE 
UNIVERSITY  OF  MICHIGAN  STATIONS 


Station/Component 

First  P 

Maximum  S 

Vertical 

Longitudinal 

Vertical 

Longitudinal 

Transverse 

Pahute  Mesa 

1.155 

1.274 

1.820 

1.470 

1.  689 

North 

0.848 

0.£;i4 

1.424 

1.355 

1.130 

East 

1.096 

1.051 

1.363 

1.213 

0.988 

Southeast 

1.862 

1.616 

1.470 

1.267 

1.514 

South 

1.197 

0.814 

1.972 

1.129 

1.585 

Mean 

Overall  Mean: 

1.191 

1.287 

1.075 

11 

1.592 

1.281 

1.352 
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Another  approach  to  the  comparison  of  the  ground  motion  generated  by  these 
two  shots  is  to  examine  the  ratio  of  the  maximum  signal  level  measured  on  the 
broadband  seismograms.  For  the  same  component  at  the  same  site  the  time  of  the 
maximum  signal  will  usually  correlate  beti^een  shots.  However,  the  maximum  sig¬ 
nal  on  one  component  will  not  necessarily  correlate  in  time  with  the  maximum 
signal  on  the  other  two  components  for  the  same  shot.  The  results  of  this 
analysis  is  sho^^n  in  Table  3.  T,  1,  and  V  are  abbreviations  for  transverse 
longitudinal  and  vertical,  respectively.  It  can  be  seen  that  the  MINE  DUST  HE 
event  was  a  factor  of  1.32  larger  than  the  DIAMDND  MINE  HE  shot  on  the  basis 
of  these  broadband  signal  ]  ivel  comparisons.  This  is  very  close  to  the  overall 
spectral  ratio  amplitude  of  1.287  shown  in  Table  2. 


TABLE  3. 

MAXIMUM  AMPLITUDE  RATIOS  MINE  DUST  HE  vs.  DIAMOND  MINE  HE 
UNIVERSITY  OF  MICHIGAN  STATIONS 


Station 

Component 

Amplitude  Ratio 

Pahute  Mesa 

L 

1.496 

(20,016  ft)* 

T 

1.413 

V 

1.585 

North 

L 

1.259 

(12,700  ft)* 

T 

1.334 

V 

1.334 

East 

L 

1.259 

(7,698  ft)* 

T 

.794 

V 

1.259 

Southeast 

L 

1.496 

(7,843  ft)* 

T 

1.413 

V 

1.334 

South 

L 

1.496 

(8,443  ft)* 

T 

1.496 

V 

1.059 

Mean  =  1.32 
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Another  measure  of  the  relative  sizes  of  the  events  is  the  energy  spectrum 
as  derived  from  theory.  The  theory  and  results  of  this  analysis  follow. 

From  Bullen  (p.  268)  we  have  the  following: 

E  =  ATT^pr^sinA/HcC^T  ^dt. 

The  symbols  and  what  each  denotes  are: 

E  total  source  energy 
p  medium  dansity 
rQ  Earth's  radius 

A  angle  between  source  and  sensor  at  the  earth's  center 
H  medium  depth 
c  wave  velocity 

C  particle  displacement  amplitude 
T  period 
t  time 

The  Integral  Is  taken  over  the  signal  of  Interest.  This  equation  describes 
the  total  energy  of  the  source  If  the  assumptions  are  made  that  the  energy  Is 
coupled  spherically  Into  the  earth  and  that  the  energy  Is  transmitted  entirely 
by  body  waves.  From  the  following; 

D  =  r  slnA,  distance  frwf  shot  to  sensor  over  the  earth's  surface 
o 

f  =  1/t,  frequency 

/V^dt  =  /(27TfC)^dt,  V  =  the  particle  velocity 
we  have 

E  =  TTpDH  /  cV^dt. 

For  our  purposes  the  Integral  Is  defined  from  the  first  arrival  time  to 
the  time  where  the  signal  returns  to  background  level.  Based  tan  an  examination 
of  the  records  the  latter  time  Is  ten  seconds  after  the  former.  We  let  t^  be 
the  time  of  the  first  motion  In  seconds,  with  t  -  0  at  the  time  of  the  shot. 
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With  this  definition  of  t.  c(t)  = 

»  \  /  ^  j 

octave  bandwidth  about  the  frequency 


and  the  energy  equation  for  a  one-third- 
f  becomes  equivalent  to 


E(f) 

irpHD^ 


to+10 

((V-:(t))Vt)dt. 


(1) 


The  actual  calculation  made  ir,  the  energy  program  is  a  constant  times 
the  digital  equivalent  of  the  right-hand  side  of  this  equation.  The  assump 
tion  is  then  made  that  for  a  given  site  the  transmitting  medium  does  not 
change  from  shot  to  shot  (i.e.,  p,  D  and  H  do  not  change). 

Since  the  ground  particle  velocity  is  sensed  in  three  mutually  orthog¬ 
onal  directions  (transverse,  longitudinal,  and  vertical),  it  may  be  expres¬ 
sed  as  the  square  root  of  the  sum  of  the  squares  of  its  components  at  any 

given  time.  Assuming  this  rule  holds  for  the  filtered  particle  velocity, 
Eq.  1  may  be  written: 


E(f) 

TTpHD^ 


to+10 


or  as  the  sum  of  three  integrals  of  the  form 
rto+lO 

4  ((V^f(t))Vt)dt 


where  V^^(t)  denotes  the  i*^^  component  of  the  particle  velocity  filtered  in 
the  one-third-octave  band  about  f  at  time  t. 

The  discrete  form  of  this  integral  is 


2  ((V,f(t  ))Vt  )  (^t).,  t  -  t+io 

joQ  J  J  jniG 


(3) 


Assuming  the  sample  interval  to  be  constant  (i.e.,  t  -t 
T)  it  may  be  removed  from  the  summation.  At  time  tj  thi 


“  ■  constant 

actual  value 
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digitized  Is 


where 


k  -  (512/3)  (10) 


-a/20 


_2 

a  “  geophone  response  IndB  re  10  cm/sec 

+  spectrum  level  correction  in  dB  (4c) 

-  total  system  gain  In  dB. 

512/3  ■  1024/6  is  the  digitizing  scale:  +3  volts  to  -3  volts  (full  record 
level  from  tape)  Is  divided  into  1024  leve’  .  (10  bits).  Note  that  k  is  in¬ 
dependent  of  time  for  a  given  frequency  and  component,  and  it  may  be  removed 
from  the  summation.  Thus  (3)  may  be  written 


T  ni  , 

T.  (U_(t,))Vc,,  t  -  t  +10. 

k^j-0  if  J  J  ra  o 


This  final  expression  is  the  digital  equivalent  of  (2)  in  terms  of  the  digitized 
data  and  the  explicitly  defined  parameters  (4). 

A  program  (Appendix  E)  was  written  to  calculate  these  values.  The  square 
roots  of  the  results  are  shown  in  Appendix  F.  The  energy  spectrum  ratios  as  a 
function  of  site  are  displayed  in  Appendix  G.  These  are  simply  point  by  point 
ratios  of  the  values  squared  given  in  the  previous  appendix.  Since  the  energy 
calculation  requires  good  data  from  all  three  components  and  since  timing  is 
critical,  those  stations  lacking  one  or  the  other  have  been  excluded.  However, 
the  data  are  sufTiclently  redundant  to  enable  quantitative  comparisons  to  be 
made  of  the  relative  energy  levels  of  the  four  shots. 
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A  number  of  qualitative  observations  may  be  made  at  this  point.  At  .'he 
North  site  the  indication  is  that  the  source  energy  peaks  at  10  Hz  and  at  the 
Pahute  Mesa  site  the  peak  has  dropped  to  4-5  Hz,  at  which  frequency  the  North 
site  contains  a  secondary  peak.  A  partial  explanation  of  this  is  that  the 
energy  calculation  does  not  take  the  absorption  of  the  medium  into  account. 

There  is  not  enough  information  available  to  estimate  an  absorption  coefficient, 
which  would  be  needed  to  account  entirely  for  the  results.  The  apparent  peaks 
also  drop  from  10  Hz  to  3-4  Hz  as  seen  by  each  station  in  turn,  looking  clock¬ 
wise  from  the  North  site  to  the  South  site.  This  again  may  be  due  to  the  lack 
of  the  inclusion  of  absorption  factors  in  the  calculation.  It  is  to  be  expected 
that  the  factors  increase  in  this  direction  because  the  geology  becomes  more 
and  more  complex  (,'aulted). 

The  Pahute  Mesa  site  data  indicate  that  the  apparent  energy  of  DIAMOND  MINE 
has  comparatively  less  high  frequency  content  than  the  other  three  shots.  This 
may  be  due  in  part  to  that  shot's  proximity  to  a  fault  along  which  the  rock 
material  (tuff)  is  likely  to  be  considerably  more  fractured.  High  frequencies 
would  be  attenuated  greater  in  this  type  of  medium.  In  short,  in  this  case 
our  assumption  that  the  medium  through  which  the  energy  travels  remains  the 
same  from  shot  to  shot  for  a  given  site  may  be  wrong.  The  azimuth  from  the 
DIAMCWD  MINE  working  point  to  the  North  site  does  not  lie  along  the  fault  as 
does  that  from  the  working  point  to  the  Pahute  Mesa  site. 

In  order  to  get  a  measure  of  the  relative  size  of  the  four  events  the 
geometric  mean  of  the  energy  ratio  values  were  calculated.  The  six  numbers 
thus  derived  are  redundant;  for  each  of  the  three  inequalities  relating  the 
sizes  of  the  events,  three  values  may  be  found.  The  three  values  were  different 
in  general,  because  all  questionable  data  were  ignored.  A  weighted  average 
(proportional  to  the  amount  of  data  used)  of  the  values  appears  in  Table  4  along 
with  other  shot  size  determinations  which  have  been  made.  These  latter  are 
means  and  ration  of  particle  velocity  measurements  and  in  order  to  compare  them 
with  the  energy  values  the  square  roots  of  the  energy  ratios  have  also  been 
included. 
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TABLE  4.  SUMMARY  OF  RELATIVE  SHOT  SIZES 


DD* 

MD* 

DD* 

DM* 

HE 

HE 

MD 

DD 

6.985 

1.233 

5.665 

6.316 

2.642- 

1.110- 

1.29 

2.38(>f 

2.513- 

1.32 

2.5 

1  A 

7  A 

mJ  • 

7  7-2 

1.7 

1.4 

2.5 

Source 


Energy  ratio 

Square  root  of  energy 
ratio 

Geometric  mean  of  spectrum 
ratios 

Geometric  mean  of  maximal 
of  broadband  signals 

Average  broadband  signalst 

Seismic  station  at  Beattvrf 

2  Hz  spectrum  ratiost 

10  Hz  spectrum  ratiost 

Coda  length  magnituict 


*DD;  DIAMOND  DUST 
MD;  MINE  DUST  HE 
DM;  DIAMOND  MINE 
HE:  DIAM(»JD  MINE  HE 

i 

See  references  3  and  4. 

The  energy  equation  which  has  been  used  is  valid  for  only  body  waves. 
Using  it  on  data  including  surface  waves  gives  results  which  are  too  large, 
since  the  surface  waves  have  a  geometric  spreading  proportional  to  the 
inverse  square  root  of  distance  whereas,  body  wave  energy  spreads  as  the 
inverse  square  of  distance.  The  seismic  motion  generated  by  piAMOND 
MINE  HE  was  predominantly  body  waves.  The  duration  of  the  surface  waves 
for  DIAMOND  DUST,  MINE  DUST  HE,  and  DIAMOND  MINE  increased  respectively, 
the  latter  having  a  coda  about  equally  divided  between  surface  waves  and 
body  waves.  Consequently,  the  energy  ratio  of  DIAMOND  MINE  to  DIAMOND  DUST, 
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for  example,  is  probably  too  large.  This  is  indicated  in  Table  4  by  the 
sign  follOTving  the  ratio.  For  the  worst  case,  DIAMOND  MINE,  the  E(f ) /  (irpD^A) 
calculation  may  be  off  as  much  as  fifteen  percent  for  the  lower  frequencies. 

It  is  recommended  that  the  energy  calculation  for  surface  waves  be  programmed 
and  used  for  those  parts  of  the  data  which  are  clearly  surface  motion  phases. 

Table  5  gives  a  breakdown  of  the  energy  ratios  derived  from  the  data  re¬ 
corded  at  the  ind'vidual  sites.  For  the  MINE  DUST  HE /DIAMOND  MINE  HE, 

DIAMOND  DUST/MINE  DUST  HE,  and  DIAMOND  DUST /DIAMOND  MINE  HE— for  those  not 
involving  DIAMOND  MINE — the  ratios  agree  among  the  sites.  This  would  indicate 
that,  even  though  the  geology  along  the  energy  paths  vary,  the  data  from  any  of 
them  give  a  good  e ’timate  of  the  energy  ratios — when  DIAMOND  MINE  is  excluded. 
This  is  what  should  be  e  q)ected  if  the  assumptions  made  in  the  above  calcula¬ 
tions  were  true. 


TABLE  5. 

ENERGY  RATIO  CALCULATIONS  BASED  ON  THE 
DATA  RECEIVED  AT  EACH  SITE 


Site 

MD* 

HE 

DM* 

MD 

DD* 

MD 

DM* 

HE 

DD* 

HE 

DM* 

DD 

Pahute  Mesa 

North 

1.463 

1.196 

5.139 

7.517 

37.93 

45.35 

East 

48.05 

Southeast 

1.587 

19.82 

4.507 

31.46 

7. 152 

4.399 

South 

1.042 

51.42 

6.516 

49.35 

6.254 

7.891 

Mean 

1.287 

36.92 

5.324 

41.29 

6.953 

5.892 

*»!:  DIAMOND  MINE 
DD:  DIAMOND  DUST 
MD:  MINE  DUST  HE 
HE:  DIAMOND  MINE 

HE 

* 
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For  those  ratios  involving  DIAMOND  MINE,  however,  the  variance  from  site 
to  site  is  quite  large,  but  consistent.  With  respect  to  the  mean  of  the  ratios, 
the  Southeast  site  data  indicates  an  average  of  49  percent  less  than  the  mean 
energy  ^rom  DIAMOND  MINE  and  the  South  site  data  indicates  31  percent  more.  The 
data  from  East  and  North  sites  Indicate  DIAMOND  MINE  to  be  a  little  larger  than 
the  means.  The  Pahute  Mesa  data  estimates  were  lower  as  was  mentioned  previously. 

One  way  to  account  for  this  azimuthal  variation  is  to  hypothesize  that  energy  was 
released  in  the  fault  next  to  which  the  DIAMOND  MINE  working  point  was  located. 

The  P  wave  radiation  pattern  generated  by  slippage  along  the  fault  nearest  the 
working  point  is  shown  in  figure  5.  The  observed  energy  ratios  agree  with  the 
hypothesis  as  can  be  seen.  The  fact  that  more  surface  wave  energy  was  generated 
during  this  event  also  supports  the  hypothesis. 

The  site  occupied  for  the  MINE  DUST  HE  event  at  a  distauca  of  225  km  did  not 
receive  adanuate  signal  levels  to  confidently  identify  the  signal  arrival.  This 
site  was  carefully  selected  for  very  low  background  levels  in  anticipation  of  the 
very  small  signal  levels.  The  average  background  level  during  the  time  of  the 
'.vent  was  approximately  5  x  10  ^  cm/sec  as  measured  over  the  data  band  from  2  Hz 
to  15  Hz.  This  level  of  "Signal  is  well  within  the  capabilities  of  the  instrumentation 
used  at  that  site.  To  overcome  this  ambient  seismic  background  would  require  a 
seismometer  array  which  could  be  processed  to  allow  the  enhancement  of  signals  from 
a  particular  direction. 


2.2.  U.  S.  GEOLXICAL  SURVEY 


Four  special  seismograph  stations  were  established  to  record  this  shot  in 
addition  to  the  regular  stations  operated  by  the  U.S.G.S.  Of  the  latter  approximately 
15  stations  were  recalibrated  and  gains  adjusted  for  this  shot.  The  locations  of 
these  stations  are  listed  in  Table  6a.  The  data  discussed  in  this  section  were 
supplied  by  Mr.  Fred  Fischer,  U.S.G.S.,  National  Center  for  Earthquake  Research, 

Menlo  Park,  California. 

The  stations  designed  India,  Tango,  Juliet,  and  Quebec  wefe  special  portable 
stations  that  were  equipped  with  magnetic  tape  recorders  and  employed  three-component 
short-period  seismometers  and  linear  arrays  of  vertical  seismometers.  The  two  most 
distant  stations  (Juliet,  A  =86  km,  and  Quebec, A  -  238  km)  received  no  detectable 
signals  from  the  shot. 
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The  U.S.G.S.  data  are  summarized  In  Tables  7  and  8.  The  velocity  model  for 
travel  time  correction  is  shown  in  Table  3b.  It  is  interesting  to  note  that  at 
station  NT18  (A  =  49.4  km,  azimuth  =  340.3  )  the  observed  ground  displacement  and 

particle  velocity  is  the  same  as  recorded  at  Beatty  (A  =  51.8  km,  azimuth  ^  254.4°). 

2.3.  NOAA/ESL 

The  Nevada  Special  Projects  Party  of  the  Earth  Sciences  Laboratories,  NatlonaJ. 
Oceanic  and  Atmospheric  Administration  operated  10  seismic  stations  to  record  this 
shot.  These  stations  were  equipped  as  shown  in  Table  9.  (Mr.  K.W.  King,  personal 
communication).  All  stations  used  at  least  one  vertical  short-period  seismometer 
and  four  stations  had  horizontal  seismometers.  Eight  of  the  10  stations  used 
magnetic  tape  recorders.  The  shot  was  well  recorded  at  three  stations,  weakly 
recorded  at  a  fourth  station,  and  not  detected  on  the  remaining  six. 

The  fact  that  the  shot  was  not  recorded  at  ETS-1,  RHN,  and  SMN  would  indicate 
that  there  were  either  strong  azimuthal  variations,  local  background  noise  masked 
the  signal,  or  the  original  record  gains  were  too  low.  As  mentioned  in  the 
previous  section,  the  signal  level  recorded  at  approximately  50  km  to  the  north 
and  to  the  west  was  equal.  As  a  distance  of  approximately  31  km  the  signal  level 
at  station  LSM  (  azimuth  191.7°)  was  a  factor  of  two  or  more  lower  than  stations 
N119  (azimuth  ==343.3°)  and  NT12  (azimuth  =  291.7°)  thus  indicating  that  there  was 
an  azimuthal  effect  with  lower  signals  recorded  to  the  south.  Comparing  stations 
CPX  and  NT16  shows  that  this  asymmetry  extends  to  the  southeast.  This  azimuthal 
asymmetry  was  also  observed  in  the  previous  shots  from  Area  16  (Ref. 3). 

2.4.  MAGNITUDE  CALCULATIONS 

The  MINE  DUST  HE  shot  did  not  generate  large  enough  signals  to  be  detected 
at  200  km  so  a  body  wave  magnitude  (nij^)  could  not  be  determined  directly.  An 
XMAG  was  computed  by  the  USGS  from  their  local  seismic  network  (see  Table  7).  This 
XMAG  is  approximately  equivalent  to  Richter's  M^^  local  magnitude.  Using  the  M^  to 
m^  conversion  discussed  in  our  earlier  report  (Ref. 3)  gives  an  m^  magnitude  of  2.2 
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TABLE  7  -  EXPLANATION 


Tht  li'sfing  of  the  so/ufton  for  the  b/ast  gufes  first  the 
pammhts  for  fk  blast  uno  -Aen  fk  o/ota.  and 

paramefers  ffora  each  senmo^raph  stafion  .  The  fol/ou)'ma 
S^niboh  are  used: 

- _  V. 

Year  month  day 
hour  (GMT)  minute  second 
LAT  N:  in  degrees  and  minutes 

LONG  W:  in  degrees  and  minutes 
DEPTH:  in  kilometers 


liAG; 

NO: 

DiM: 

GAP : 

RMS; 

ERH: 

ERZ: 

Q: 


duration  magnitude,  approximately  equal  to  Richter 
number  of  first-arrival  times  actually  used  in  the  location, 
epicentral  distance,  in  kilometers,  to  the  nearest  station, 
largest  azimuthal  separation,  in  degrees,  between  stations, 
root  mean  square  error  of  time  residuals  in  seconds, 
standard  error  of  the  epicenter  in  kilometers, 
standard  error  of  the  focal  depth,  in  kilometers. 

solution  quality  of  the  hypocenter;  A  =  excellent,  B  =  good, 
C  =  fair,  D  =  poor. 


For  Jie  individual  station  parameters  the  following  symbols  are  used 
•  station  code  name 

D^ST:  epicentral  distance  to  the  station,  in  kilometers. 

AZM:  azimuthal  angle  to  the  station,  from  0°  north. 

AIN:  angle  of  inclination  at  which  the  first-arrival  ray  leaves 

the  earthquake  focus. 


I 


Table  7  Explanation  (cont) 


PRMK: 

HRMN: 

P-SEC: 

TP03S: 

TPCAL; 

DELAY: 

P-RES : 

P-WT: 


I  =  impulsive  first  arrival;  E  =  emergent;  P  =  phase  remark; 
U  =  compressional  first  motion;  D  =  dilatational  first 
motion. 

hour  minute  (GMT) . 
first-arrival  time, 
observed  P-wave  travel  time, 
calculated  P-wavc  travel  time. 

station  delay,  introduced  at  the  beginning  of  the  program 
to  compensate  for  local  variations  in  the  velocity  model, 
travel  time  residual,  in  seconds;  (TPOBS)- (TPCAL)- (DELAY)  = 
P-RES. 

weight  assigned  to  station  reading  in  the  location;  a 
decreasing  function  of  distance. 


25 


MJ 


5c  “ 


> 

V) 

1 

•o 

V  ^ 

'o 

5 

'o 

o  2 

Ve 

k 

X 

X  X 

C> 

c: 

>0 

fO 

ro 

:^ 

N 

'o  'o  'c 


>  tn  Q 

r>J  ®i  -C  < 


^  •* 
Q  O 


«r 

<0 


» 

X 

^  o 

11 

k 

|( 

•M 

It 

'V.  -o 

§  s: 

:i 

1 1 

Ki 

o 

11 

^  s 

1 

w 

i-t 

00 

!) 

Vl 

•ft 

P: 

'K 

Of 

o 

^  'Ji! 

v^  vn 

o 

§ 

S 

1  ’ 

5 

i\  II 

5  u.  o 
>  "5:  ^ 


>  >  \ 
o  S  '- 

'•  K 

T-l  >  '<■ 

:r 

U.  ^  ^ 

5  Vf)  vS 


O  O'  X 
rr  fy 

m  S*  a. 

ft. 


t- 

«^'5, 

m  :a  « 

s  5  < 


e? 

UJ 

KJ  ^ 

I  u<l 

^  5 

P  to 


w  '* 

>  00 


Q  O  ^ 

N  w  CM 


!o  ^  lo  o 

cS  ^  ci  $  cJ 


**  ^  ^  o 

fft  ir 


o  ^.  N  ^ 

s 

P>  fo  ^  ^ 


'cs  -« 
H  3 


o  2  ^  o 

^  ri 

^  .  .  . 

VU  ^  o  o 

Q. 


A  2 

o 

«  i" 

A  ft<  V- 
R  U1  ? 
R  v>  V- 
lo  ^  M) 


5  ^  -2 


to  00 

v:  k 


vo^tu  vuowJDif  »/cftw($  bvtiunS9> 


(NVIRONMINTAl 

RESEARCH  INSTITUTE 
OF  MICHIGAN 


FORMERLY  WILLOW  RUN  LABORATORIES.  THE  UNIVERSITY  OF  MICHIGAN 


P  0  BOX  618  •  ANN  ARBOR  •  MICHIGAN  •  48107 


PHONE  (313)  483  0500 


Which  is  obviously  too  high.  Table  10  summarizes  the  coda  length  magnitudes 
(Mp  computed  from  the  USGS  network  (Mr.  Fred  Fischer,  personal  communication)  for 
the  four  events  of  this  series.  The  MINE  DUST  HE  shot  is  0.23  magnitude  units 
larger  than  the  earlier  HE  shot  and  0.15  magnitude  units  lower  than  DIAMOND  DUST. 
These  M^  values  are  believed  to  be  significant  insofar  as  their  relationship  to 
each  other.  However,  the  dilemma  of  computing  an  accurate  body  wave  magnitude  for 
these  small  events  has  still  not  been  resolved.  The  body  wave  magnitude  is 
important  since  that  is  the  basis  on  which  seismologlcal  decoupling  estimates  have 

been  made. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

This  FOR  complements  the  previously  Issuer  technical  report  (ref.  5)  for 
the  MINE  DUST  HE  event  and  contains  the  final  results  of  the  analyses  which 
have  been  performed  on  all  of  the  MIGHTY  MITE  events  but  which  have  not  been 
reported  previously  (refs.  3,  4,  5).  It  was  originally  planned  to  Instrument 
the  same  sites  In  the  same  way  for  the  MINE  DUST  nuclear  event.  Construction 
difficulties  have  postponed  this  event  Indefinitely.  The  following  conclusions 
are  a  summary  of  the  results  reported  herein. 

1.  The  geometric  mean  of  the  spectral  ratios  for  the  first  compresslonal 
wave  arrival  and  the  maximum  shear  wave  re  :orded  at  the  five  close—ln  stations 
Indicate  that  the  MINE  DUST  HE  shot  was  a  factor  of  1.287  larger  than  the 
DIAMOND  MINE  HE  shot. 

2.  The  ground  velocity  amplitudes  measured  from  the  broadband  seismograms 
Indicate  that  the  MINE  DUST  HE  shot  was,  at  the  mean,  a  factor  of  1.32  larger 
than  the  DIAMOND  MINE  HE  shot. 

3.  Coda  length  magnitudes  calculated  for  MINE  DUST  HE  and  DIAMOND  MINE  HE 
Indicate  that  the  MINE  DUST  HE  shot  was  0.23  magnitude  units  larger. 

4.  The  MINE  DUST  HE  shot  did  not  generate  large  enough  signals  to  deter¬ 
mine  body  wave  magnitude  directly. 

5.  The  geometric  mean  of  an  energy  ratio  measure  Indicated  that  MINE 
DUST  HE  was  1.23  times  larger  than  DIAMOND  MINE  HE;  that  DIAMOND  DUST  was 
5.67  times  larger  than  MINE  DUST  HE;  that  DIAMOND  MINE  was  6.32  times  lasger 
than  DIAMOND  DUST.  These  values  are  energy  ratio  values,  so  they  are  com¬ 
parable  to  the  squares  of  the  velocity  ratio  values  given  above. 

6.  There  Is  evidence  that  tectonic  strain  release  may  have  occurred  at 
the  time  of  DIAMOND  MINE. 

The  following  recommendations  are  made  based  on  these  conditions. 

1.  Body  wave  magnitude  measurements  have  not  been  possible  because  suffi¬ 
cient  data  have  not  been  gathered  at  a  distance  of  200  km  or  more.  An 
alternative  to  having  a  shot  In  the  same  area  which  would  be  large  enough  to 
record  at  that  distance— this  not  being  feasible— Is  having  a  smaller  shot 
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and  deploying  a  sensor  array  at  the  appropriate  distance.  Summing  the  sensor 
outputs  or  other  array  processing  techniques  could  then  be  employed  to  dis¬ 
tinguish  the  signal  of  Interest  from  the  background. 

2.  An  analytic  approach  to  measuring  the  body  wave  magnitude  differences 
between  the  shots  may  be  possible  by  extrapolating  on  the  data  already  re¬ 
corded.  This  will  require:  a)  the  preshaping  of  the  P-wave  train  data  with  a 
filter  matched  to  the  response  of  the  World-Wide  Standard  Seismograph  Network 
short-period  Instruments,  and  b)  computing  the  square  root  of  the  P-wave 
energy  for  each  component  taken  separately  and  in  combination.  The  results 
can  then  be  compared  with  existing  data  on  the  P^  branch  of  the  travel  time 
curve. 

3.  As  mentioned  in  Section  2.1,  the  energy  ratio  calculations  can  be 
made  more  accurate  by  using  a  separate  equation  for  the  surface-wave  phases 
of  the  coda.  This  will  require:  a)  writing  another  program  comparable  to  the 
one  In  Appendix  E  for  the  body  wave  energy  calculations,  b)  obtaining  para¬ 
meters  from  the  records  which  would  define  a  partition  of  the  data  into  surface 
and  body  phases,  and  c)  rerunning  the  programs  using  the  data  which  has  already 
been  digitized. 

4.  Tectonic  strain  releiise  apparently  occurred  at  the  time  of  DIAMOND 
MINE,  which  further  confuses  the  determination  of  its  relative  energy.  If 
there  were  strain  release  the  total  energy  radiated  would  be  the  sum  of  the 
energy  from  both  sources  less  the  amount  of  energy  from  DIAMOND  MINE  which 
was  transformed  to  the  mechanical  work  needed  to  trigger  the  strain  release. 

The  data  and  theory  are  sufficient  to  get  a  first  order  estimate  of  the  pro¬ 
portion  of  the  energy  due  to  tectonic  strain  release.  This  could  be  accomplished 
by  firing  the  observed  data  to  the  theoretically  derived  radiation  patterns 
generated  by  the  mechanism  hypothesized  in  Section  2.1. 
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APPENDIX  A 

MINE  DUST  HE  Spectral  Analyses 
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Figure  A-4.  MINE  DUST  HE,  NTS,  East  Site  (Pole  Line  Road) 

(Longitudinal) 
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Figure  A-5.  MINE  DUST  HE,  NTS,  East  Site  (Pole  Line  Road) 

(Transverse) 


38 


P  max . 

S  max  — — 

Avr.  Background  - 


Figure  A-6.  MINE  DUST  HE,  NTS,  East  Site  (Pole  Line  Road) 

(Vertical) 
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Figure  A-9.  MINE  DUST  HE,  NTS,  Southeast  Site  (Vertical) 
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Figure  A-10.  MINE  DUST  HE,  NTS,  South  Site  (Longitudinal) 
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Figure  A-11. 


MINE  DUST  HE,  NTS,  South  Site  (Transverse) 
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Figure  A-12.  MINE  DUST  HE,  NTS,  South  Site  (Vertical) 
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Figure  B-9.  South  Site,  NTS,  First  P  (Vertical) 


CNVIRONMCNTAL 

RESEARCH  INSTITUTE 
iiJ  OF  MICHIGAN 


FORMERLY  WILLOW  RUN  LABORATORIES.  THE  UNIVERSITY  OF  MICHIGAN 


P  0  BOX  618  •  ANN  ARBOR  •  MICHIGAN  •  48107 


PHONE  (313)  483-0500 


APPENDIX  C 

Spectral  Ratios  of  MINE  DUST  HE/DIAMOND  MINE  HE 
as  a  Function  of  Site  and  Wave  Type 


Figure  C-4.  MINE  DUST  HE/DIAMOND 
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Figure  C-9.  MINE  DUST  HE/DIAMOND  MINE  HE,  S  ,  Vertical 
®  *  max’ 
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APPENDIX  D 

Spectral  Ratios  of  DIAMOND  MINE/DtaMOND  DUST 
as  a  Function  of  Site  and  Wave  Type 
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Figure  D-1.  DIAMOND  MINE/DIAMOND  DUST,  First  P,  Longitudinal 


Figure  D-2.  DIAMOND  MINE/DIAMOND  DUST,  First  P,  Transverse 


Figure  D-3.  DIAMOND  MINE/DIAMOND  DUST,  First  P,  Vertical 
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Figure  D-5.  DIAMOND  MINE/DIAMOND  DUST 


Figure  D-6.  DIAMOND  MINE/DIAMOND  DUST,  Vertical 
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Figure  D-7,  DIAMOND  MINE/DIAMOND  DUST,  S  ,  Longitudinal 
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FORTRAN-SABR  Program  for  Body  Wave 
Energy  Calculations 


W5PIL 

m\  C . PROGRAM  MAME!  E3718*  DATE  6/19/72 

002  C . PARTICLE  ESERGY  PRO(3?AM. 

003  C . USES  DATA  TAPE  #2*  R  CH.  MULTIPLEXED. 

004  COMMOM  IDRIF*IBLIF>IDX 

005  DIMENSIOM  IDBUEC5)»IBUE(1024#2)#IDXC9)*AC3*8)»T<8)»T0<8>» 

006  1EC8)*EEC8)*SHIFTC8) 

007  CALL  lOPENC ’PTRSD) 

008  S  $TA  /LOAD  -1 

009  PAUSE 

010  S  CLA  /CLAER  THE  ACC. 

011  1  DO  3  1=1*3 

012  3  READ(4*100)IDXCI)*CA(I*J)*J=1*8) 

013  RFAn(4*101)T0*TS*WW 

014  DO  4  1=1,8 

015  A  FECI)=0.0 

016  DELTA=8./5000. 

017  C . START  LOOP  FOR  THREE  COMPONENTS:  VERT.  LONG.  «  TRANS. 

018  DO  70  IDN=1,3 

019  WRITFC1,102)IDXCIDN) 

020  READ(4,104)SHIFT 

021  C . ZERO  EC8)  AND  INITIALIZE  T<J)S 

02?>  ISKIP=(TS-T0(8))’*'625./256. 

023  II=ISKIP+2 

024  DIF=FLOATaSKIP)*DELTA 

025  DO  5  1=1,8 

026  T<n=T0CI)+DIF 

027  5  E(I)=0.0 

0?8  CALL  BIIFIL(TDN,II) 

029  C . READ  IN  FIRST  RECORD  NEEDED 

5,3Pj  . START  TRANSFORMING  OF  IBUF 

031  10  DO  50  IS1IB=1,2 

032  DO  50  1=1,128 

033  DO  25  J=l,8 

034  C . CHECK  IF  OUTSIDE  WINDOW  WIDTH,  IF  SO  C=0.0 

035  TT=T(J) +.00005 

036  IF(TT-TS)25,14,14 

037  14  IF(TT-TS-WW)1 7,1 7,25 

038  C . FORM  PARTIAL  SUMMATION  FOR  EACH  OF  8  PARTS 

039  17  II=8*(I-1)+J 

040  X=IBUF(II,ISUR) 

041  X=X-$HIFT<J) 

042  E(J>=E<J)+(X*X)/T(J) 

043  25  CONTINUE 

044  TT=T<1 )♦. 00005 

045  IFCTT-TS-WW)40, 40,30 

046  C . FORM  PARTIAL  SUM  OF  EE 

047  30  DO  35  L=l,8 

048  X=ACIDN,L> 

049  ECL>=E(L>*X*X 

050  35  EECL)=EECL)+E(L) 
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1 


50* /L 

p,50  35  FE(L)=EFCL)+ECL) 

051  WRITEC1»103)E 

052  GO  TO  70 

053  40  DO  45  L=1»8 

. . ipoate  time  COUMTERS 

055  45  TCL)=TCL)+DELTA 

056  50  COVTI^JUE 

057  C . REFILL  IRUF  WITH  NEXT  RECORD 

055  CALL  BtJFIL(0*ISKIP) 

059  GO  TO  10 

0«!J  70  COMTI>>UE 

061  V,'RITE(1*103)EE 

062  C . .START  ALL  OVER  AGAIN 

0fi3  C......ANn  rx3  A>JOTHER  ID 

0A4  GO  TO  1 

065  100  FORMAT(I4.8F10.8) 

066  101  FORMAT(9E12.8*F5.0> 

067  102  FORMAT(/*ID  #  *.14/) 

068  103  F0RMAT<4E16.8/) 

069  104  E0RMAT<8E16.8) 

070  END 
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Energy  Program  Results  as  a  Function 
of  Site  and  Event 
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Figure  F-4 
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APPENDIX  G 

Energy  Spectrum  Ratios 
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Figure  G-3.  DIAMOND  MINE/DIAMOND  MINE  HE 
spectrum  ratios 
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Figure  G-A.  DIAMOND  MINE/DIAMOND  DUST  energy 
spectrum  ratios 


